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ABSTRACT 

We present the first systematic study of X-ray flare candidates in short gamma-ray 
bursts (SGRBs) exploiting the large 6-year Swift database with the aim to constrain 
the physical nature of such fluctuations. We find that flare candidates appear in differ- 
ent types of SGRB host galaxy environments and show no clear correlation with the 
X-ray afterglow lifetime; flare candidates are detected both in SGRBs with a bright 
extended emission in the soft 7-rays and in SGRBs which do not show such compo- 
nent. We furthermore show that SGRB X-ray flare candidates only partially share the 
set of observational properties of long GRB (LGRB) flares. In particular, the main 
parameter driving the duration evolution of X-ray variability episodes in both classes 
is found to be the elapsed time from the explosion, with very limited dependence on 
the different progenitors, environments, central engine life-times, prompt variability 
time-scales and energy budgets. On the contrary, SGRB flare candidates significantly 
differ from LGRB flares in terms of peak luminosity, isotropic energy, flare-to-prompt 
luminosity ratio and relative variability flux. However, these differences disappear 
when the central engine time-scales and energy budget are accounted for, suggesting 
that (i) flare candidates and prompt pulses in SGRBs likely have a common origin; 
(ii) similar dissipation and/or emission mechanisms are responsible for the prompt 
and flare emission in long and short GRBs, with SGRBs being less energetic albeit 
faster evolving versions of the long class. Finally, we show that in strict analogy to the 
SGRB prompt emission, flares candidates fall off the lag-luminosity relation defined 
by LGRBs, thus strengthening the SGRB flare-prompt pulse connection. 
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1 INTRODUCTION 

With an isotropic peak luminosity up to 10 54 ergs -1 , 
gamma-ray bursts (GRBs) are the brightest objects in the 
7-ray sky during their short lives (At ~ 0.1 — 100 s). 
Their duiation-spectral hardness distribution gives evidence 
for the presence of two classes (Kouveliotou et al. 1993): 
long and short GRBs (observed duration longer and shorter 
than 2 s, respectively), with short bursts appearing slightly 
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harder. The dichotomy in the duration-hardness dimensions 
suggested separate progenitor populations. However, until 
a few years ago, the distances, energy and environments of 
SGRBs (short GRBs) remained highly uncertain due to the 
poor localisation. 

The breakthrough in the study of SGRBs occurred 
thanks to the rapid slew capabilities of the Swift spacecraft 
(Gehrels et al. 2004) which allowed spectroscopic observe 
tions to be performed at very early times. These observer 
tions revealed that SGRBs are cosmological, with prompt 
luminosities comparable to LGRBs albeit significantly less 
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energetic; with similar afterglows (Nysewander et al 2009) 
but residing in completely different environments. In sharp 
contrast to LGRBs, short bursts have been localised both in 
early- type and late- type host galaxies (see Berger 2011 and 
references therein), pointing to an old progenitor popula- 
tion. The detection of supemovae associated to LGRBs (see 
e.g. Kulkarni et al. 1998; Stanek et al. 2003; Fruchter et al. 
2006 and references therein) provided instead support to 
models invoking young stellar progenitors. According to 
the standard scenario, LGRBs originate from the collapse 
of rapidly-rotating, massive stars (MacFadyen & Woosley 
1999), while SGRBs are believed to result from the coales- 
cence of a binary system of compact objects (neutron star 
plus neutron star NS+NS or neutron star plus black hole 
NS+BH, Paczynski 1986; Eichler et al. 1989; Narayan et al. 
1992). 

Despite fundamental theoretical and observational 
progress, the nature of SGRB progenitors remains elu- 
sive. Numerical simulations show that the active stage of 
a NS+NS merger typically lasts ~ (0.01 — 0.1) s (see e.g. 
Nakar et al. 2007 and references therein 1 ): material ejected 
during the merger is expected to accrete on time-scales of 
the order of 1 — 10 s (the exact value depending on the ac- 
creting disk viscosity parameter and details of the ejection 
process). Thus, the detection of central engine activity on 
time-scales much longer than the usual dynamical or even 
viscous time-scales would challenge the currently accepted 
scenario (see Nakar et al. 2007 for a recent review). 

Long-lasting (At 10 s), soft energy tails detected in 
several SGRBs during their prompt emission (the so-called 
extended emission, see Norris et al. 2010a and references 
therein) represent such a case and pose severe constraints 
to existing models, especially when energetically dominat- 
ing with respect to the primary burst (Perley 2009). The 
same is true for the recently discovered presence of precur- 
sors (Troja et al. 2010). Equally challenging would be the 
detection of late-time central engine activity in the form of 
flares superimposed over the smooth SGRB X-ray afterglow. 

Flares are currently detected in ~ 30% of long 
GRBs X-ray afterglows (Chincarini et al. 2010) as fast- 
rise exponential-decay features whose spectral and temporal 
properties have been demonstrated to show a strict analogy 
to LGRB prompt pulses (Margutti et al. 2010b): this find- 
ing suggested that flares might originate from re-activations 
of the LGRB central engine. Several ideas on how to explain 
the possible presence of flares in short GRBs have been ex- 
plored as well: the fragmentation of the outer parts of an 
hyper-accreting disk around the newly formed black hole as 
a result of gravitational instabilities could potentially lead 
to large-amplitude variations of the central engine output 
of both long and short GRBs (Perna et al. 2006). Alterna- 
tively, the late-time accretion of material launched into ec- 
centric but gravitationally bound orbits during the compact 
binary merger could provide the fuel to revive the central 
engine activity (Rosswog 2007). The long term evolution of 
debris following the tidal disruption of compact objects has 
been identified by Lee et al. (2009) as a feasible mechanism 
to produce flares. Finally, as an alternative in the context of 
accretion-powered models, magnetic halting may also give 


1 In a recent study Rezzolla et al. (2011) found At ~ 0.3 s. 


rise to secondary episodes of delayed activity as suggested 
by Proga & Zhang (2006). However, the observational prop- 
erties of flares in SGRBs have not been determined, yet, so 
that it is at the moment unclear if any of these models would 
be able to explain the observations. 

While SGRB X-ray light-curves clearly show temporal 
variability superimposed over a smooth decay, the presence 
of real flares in short bursts is questionable. In particular, it 
is at the moment unclear if what is currently identified as 
short GRB flare emission (see e.g. La Parola et al. 2006 for 
GRB 051210) quantitatively shares the very same properties 
of the population of long GRB flares: are there fast varying 
At/t -C 1, prominent temporal features in the afterglow of 
SGRBs with properties reminiscent of the long GRB flaring 
emission? Do SGRB flare candidates follow the entire set of 
relations found from the analysis of real flares in long bursts? 
In particular: is the evolution of their temporal and energetic 
properties compatible with the flare-like behaviour identified 
by Chincarini et al. (2010)? What is the typical amount of 
energy released during such episodes of variability? Is there 
any link between the late-time variability which appears in 
the X-ray afterglow of SGRBs and their prompt emission? 
Negligible spectral lag is a defining characteristic of SGRB 
prompt pulses: is this picture still valid when considering 
their late-time variability? 

Prompted by this set of still open questions, we present 
the first systematic study of X-ray flare candidates in short 
GRBs, taking advantage from the large Swift 6-year datar 
base. Through a homogeneous temporal and spectral anal- 
ysis of the widest sample of SGRB light-curves available at 
the time of writing, this study allows us to perform a one-to- 
one comparison with the properties of X-ray flares detected 
in long bursts (Chincarini et al. 2010, Margutti et al. 2010b, 
Bernardini et al. 2011, Margutti et al. 2011). The primary 
goal of this paper is to observationally constrain the origin 
of SGRB flare can didates providing the reader with a com- 
plete picture of their properties. 

This work is organised as follows: the sample selection 
and data reduction is presented in Section 2. Results are de- 
scribed in Section 3 and discussed in Section 4. Conclusions 
are drawn in Section 5. 

The GRB phenomenology is presented in the observer 
frame. Isotropic equivalent luminosities and energies are 
listed. The observer frame 0.3-10 keV energy band is adopted 
unless specified. The zero time is assumed to be the trig- 
ger time. We use the notation: V^ GRB (Y^ GRB ) to indi- 
cate that Y refers to the flare (prompt) emission of SGRBs 
(LGRBs). All the quoted uncertainties are given at 68% con- 
fidence level (c.l.). Standard cosmological quantities have 
been adopted: Ho = 70Kms -1 Mpc -1 , = 0.7, Om — 0.3. 


2 SAMPLE SELECTION AND DATA 
ANALYSIS 

We select the short GRBs detected by the Swift Burst Alert 
Telescope (BAT, BartheLrny et al. 2005a) and promptly re- 
pointed by the Swift X-Ray Telescope (XRT, Burrows et al. 
2005) between April 2005 and February 2011. The short na- 
ture of each event is established using the combined infor- 
mation from the duration, hardness and spectral lag of its 
prompt 7-ray emission: a prompt 7-ray duration Too <2s 
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coupled to a hard 7-ray emission with photon index T < 1.5 
and a negligible 7-ray spectral lag are considered in- 
dicative of a short GRB nature (see Thble 1). The mor- 
phology of the host galaxy is also used as an additional 
indicator, when available. The final sample comprises 60 
SGRBs. The presence of X-ray variability in each SGRB is 
investigated following the method by Margutti et al. (2011), 
used to determine the presence of flares in long GRBs. Only 
GRBs showing fluctuations with a minimum 2 a 2 signif- 
icance with respect to the continuum have been consid- 
ered in the following analysis. This procedure automati- 
cally identifies the best time intervals to be searched for 
the presence of X-ray flare candidates in SGRBs. Out 
of ~60 Swift SGRBs, 8 satisfy the variability require- 
ment above (Table l) 3 . Notably, the sample includes the 
unique 2 SGRBs with secure early-type host identification: 
GRB 050724 (Barthelmy et al. 2005b) and GRB100117A 
(Fong et al. 2010). In three cases (GRB 050724, GRB 070724 
and GRB 071227, in boldface in Table 1 ) an extended emis- 
sion (EE) has been detected in the soft gammarray en- 
ergy range after the short hard spike (Norris et al. 2010a; 
Norris et al. 2011). In the other cases, an upper limit on the 
EE to IPG (Initial Pulse Complex) intensity ratio (Ri a t = 
EEint/IPCfat) has been provided by Norris et al. (2010a): 
for the sample of events without EE the upper limit on 
Rint is found to be a factor > 10 below the typical Ri nt 
of SGRBs with detected EE (Table 1, column 7). Finally, 
GRB 1008 16 A has not been included in the sample in spite 
of its Tgo = 2.9 ±0.6 s (Markwardt et al. 2010) since the low 
statistics prevents the 7-ray lag analysis from giving defini- 
tive results on its possible short nature (Norris et al, 2010b). 
The burst is however considered a SGRB in Norris et al. 
( 2011 ). 


2.1 Swift-BAT data analysis 

BAT data have been processed using standard Swift-BAT 
analysis tools within heasoft (v. 6.10). In particular, the 
BATGRBPRODUCT script has been used to generate event lists 
and quality maps necessary to construct 4 ms mask- weighted 
and background-subtracted light-curves in the 50-100 keV 
and 100-200 keV anergy bands. The ground-refined coor- 
dinates provided by the BAT-refined circulars have been 
adopted; standard filtering and screening criteria have been 
applied. 


2.2 Swift-XRT data analysis 

XRT data have been processed with the latest heasoft re- 
lease available at the time of writing (v. 6.10) and corre- 
sponding calibration files: standard filtering and screening 

2 A 3cr threshold would only exclude GRB 051210, where the fluc- 
tuation has a significance of ~ 2.8cr. 

3 The percentage of SGRBs with variable XRT light-curve 
8/60 ~ 13% is much less than the ~ 30% of LGRBs showing 
flares (Chincarini et al. 2010). This result suggests that the per- 
centage of SGRB light-curves with variability superimposed is 
lower than in LGRBs. However, the lower statistics characteris- 
ing the SGRB curves prevents us from drawing firm conclusions. 
This topic will be addressed in a separate work. 



Figure 1. Up-per panel: 0.3-10 keV count-rate light-curve of 
GRB 1001 17A. Black solid line: continuous X-ray emission un- 
derlying the flare candidates computed as described in Section 
2.2; dashed lines: best- fitting flare candidate emission; red solid 
line: best estimate of the total emission. The vertical dot-dashed 
lines mark the flare candidate onset times. Inset: Complete Swift- 
XRT light-curve. The yellow filled area marks the time window for 
the computation of the CCF lag (Sect. 2.3). Middle panel: hard- 
ness ratio (HR) evolution with time; the HR is computed between 
1.5-10 keV (hard band) and 0.3-1. 5 keV (soft band). Lower panel: 
Spectral photon index evolution with time as calculated by Evans 
et al., 2010. 

criteria have been applied. Pile-up corrections have been ap- 
plied when necessary (Romano et al. 2006; Vaughan et al. 
2006). Count-rate light-curves have been extracted in the 
total XRT 0.3-10 keV energy band as well as in the 0.3-1 
keV, 3-10 keV, 0.3-1.5 keV, 1.5-10 keV and 4-10 keV en- 
ergy bands. The 0,3-10 keV count-rate light-curves have 
been re-binned at a minimum signal-to-noise ratio SN~4 and 
then searched for statistically significant temporal variabil- 
ity superimposed over a smooth afterglow decay. A two-step 
procedure has been followed: first the smooth continuum 
contribution has been determined applying the method by 
Margutti et al. (2011). A simple power-law or a smoothly 
joined broken power-law model is adopted (black solid line 
of Fig. 1). As a second step, the properties of statistically 
significant fluctuations with respect to the continuum have 
been determined adding a number of Norris et al. (2005) 
profiles to the best fitting continuum model. The best fitting 
Norris et al. (2005) profiles constitute the sample of X-ray 
flare candidates of SGRBs analysed in this work. Figure 1 
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shows GRB100117A as an example: 3 distinct episodes of 
variability have been identified. The best fitting parameters 
of the entire sample are listed in Table Al. The choice of the 
Norris et al. (2005) profile allows us to perform a one-to-one 
comparison with the properties of X-ray flares and prompt 
pulses in LGRBs (Chincarini et ah 2010; Bernardini et al. 
2011): Fig. 2 shows the evolution of the SGRB flare candi- 
dates width with time compared to LGRB flares. 

The evolution of the spectral properties of any source 
can be constrained through the analysis of its hardness ratio 
(HR), which is here defined as HR = i° 5 tv) ■ A 

different binning with respect to the total 0.3-10 keV light- 
curve has been used for the 1.5-10 keV and 0.3- 1.5 keV 
light-curves to improve the HR signal- to- noise ratio. The 
temporal evolution of the spectral photon index T has been 
calculated by Evans 2010. Results are portrayed in Fig. 1, 
Bl, B2, B3 and B4. 

Count-rate light-curves have been converted into flux 
and luminosity (when possible) curves using the spectral 
information derived from a time-resolved spectral analysis 
where the spectral evolution of the source, if present, is prop- 
erly accounted for. This procedure allows us to convert the 
best fitting peak count rates of the X-ray flare candidates 
{A parameter of Table Al) into 0.3-10 keV peak luminosi- 
ties Tpk,F B when the redshift of the source is known. Figure 
3 shows the evolution of the SGRB flare candidate 
with time compared to the results obtained for LGRB flares 
by Chincarini et al. (2010); a comparison of the two distri- 
butions can be found in Fig. 4. The isotropic equivalent en- 
ergy JSjgo p B has been determined integrating the best-fitting 
Norris et al. (2005) luminosity-calibrated profiles from the 
onset time (t s ) to t s + lOOw (where w is the flare candidate 
width) . The uncertainty arising from the spectral calibration 
has been propagated following standard practice into the fi- 
nal Igg? 3 and £?iso t F B uncertainties listed in Tables 1 and 
Al. We refer the reader to Margutti et al. (2010a) for details 
on the light-curves and spectra extraction. Figure 5 shows 
the SGRB flare candidates .Ejlo.iP distribution compared 
to the values determined for LGRB flares, as computed by 
Chincarini et al. (2010). 


2.3 Spectral time-lag computation 

The spectral lag is the time difference between the arrival 
of high-energy and low-energy photons. For each GRB, the 
X-ray and 7-ray spectral lags (r^ and T ilg’ respectively) 
and associated uncertainties have been determined using a 
cross-correlation function (CCF) analysis. The CCF analysis 
requires the observations to have a fractional exposure equal 
to 1: this requirement excludes most of XRT observations 
taken in photon counting (PC) mode 4 . Among these, the 
late- time (t ~ 5 x 10 4 s) re-brightening of GRB 050724 (Fig. 
Bl). We closely follow the prescriptions by Stamatikos et al. 


4 Swift-XKT automatically switches to the PC observing mode 
for count-rates below a few count s -1 to minimise the presence of 
pile-up. In PC mode, it is not uncommon to have short time inter- 
vals of no observation even during a single orbit. While the light- 
curve and spectra extraction procedures are basically insensitive 
to these short pauses, the CCF analysis would give un-reliable 
results. 



Figure 2. Rest frame width vs. peak time relation for LGRB 
early- time flares from Chincarini et al., 2010 (open triangles) 
and SGRB flare candidates with and without extended emis- 
sion (red open find filled circles, respectively) . Inset: Complete 
view of the w/(l + z) vs, tpk/(l + z) relation established by 
LGRB flares obtained joining the data from Chincarini et al., 
2010 and Bernardini et al,, 2011. The blue dashed line in both 
plots marks the best fitting relation calculated on LGRB flares: 

\ where w and t p k are mea- 
sured in seconds. 



Rest frame t pk (s) 

Figure 3. 0.3-10 keV peak luminosity evolution with time 
for LGRB flares (black open triangles, from Chincarini et al. 
2010) and SGRB flare candidates with and without extended 
emission (red open and filled circles, respectively). Blue dashed 
line: best fitting power-law model for LGRB flares: = 

/ \ -2.7±0.5 

10 54.8±o.4 f _£k J and extrinsic scatter a = 0.73 ± 0.08. 

Orange dot-dashed line: best fitting decay re-normalised 

by a factor 100 to match the observed SGRB flare candidates 

rSGRB 

■kpk,F ’ 

(2009) and Ukwatta et al. (2010) for the CCF computation: 
in particular, each CCF peak has been fitted using a third 
order polynomial; the number of points to be fitted around 
the CCF peak has been allowed to vary from case to case 
with the possibility to specify asymmetric intervals around 
the peak. In our analysis, a positive spectral lag is obtained 
if high energy photons lead low energy photons. 

The lag extraction is sensitive to a number of parame- 
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Figure 4. 0.3-10 keV (observer frame) isotropic equivalent peak 
luminosity Lp£p B of LGRB flares from Chincarini et al., 2010 
(hatched histogram) compared to SGRB flare candidates (filled 
histogram); two vertical arrows mark the position of flare can- 
didates in SGRBs with extended emission. The vertical dashed 
lines mark the median values of the two distributions: L^G p B ^ 
10 50 1 erg3 _1 ; Lp kt p B ~ 10 478 ergs -1 . Inset: 1-10000 keV rest 
frame isotropic equivalent peak luminosity distribution of the long 
(Nava et al., 2008) and short GRBs prompt emission (Ghirlanda 
et al., 200); 2010), with median value: Lp^p ~ 10 52,3 ergs -1 . 



Lo «[ E too.r/erg] 

Figure 5. 0.3-10 keV (observer frame) isotropic equivalent energy 
of LGRB flares from Chincarini et al., 2010 (hatched his- 
togram) compared to SGRB flare candidates (filled histogram); 
two vertical arrows mark the position of flare candidates in 
SGRBs with extended emission. Vertical dashed lines: median 
Eg™ ~ 10 51 erg and ~ 10 49 erg values. Inset: 1-10000 

keV rest frame isotropic equivalent energy distribution for the 
prompt emission of the widest samples of long (hatched his- 
togram) and short GRBs (filled histogram) with firm spectral pa- 
rameter estimates at the time of writing (Amati et al., in prep.), 
with median values: ~ 10 53 1 erg; ~ 10 51 * 5 erg. 


ters: energy band pass of each comparative light-curve, tem- 
poral bin resolution, signal-to-noise ratio and presence of 
background emission (i.e. in X-rays, the smooth X-ray decay 
underlying the time- variable signal). Fbr the prompt 7- ray 
phase, the lag has been calculated using 4 ms light-curves 
(Section 2.1) in the 50-100 keV and 100-200 keV energy 


bands. Time intervals covered by extended emission have 
been excluded. This allows us to perform a direct compari- 
son with the time-lag values obtained for LGRBs observed 
by BAT (Ukwatta et al. 2010). Results are listed in Table 1, 
column 5: the 8 SGRBs exhibit negligible ^ag- 
in the X-rays the situation is complicated by the pres- 
ence of a smoothly declining afterglow emission underly- 
ing the episodes of possible activity (see e.g. black solid 
line of Fig. 1). Choices of re-binning time-scales, energy 
bands and/or temporal intervals giving origin to correla- 
tion values (CCF peak) < 0.4 have been discarded. The 
choice of the energy bands to be compared is limited by 
the XRT 0.3 — lOkeV coverage. For each SGRB, the X-ray 
time lag T lag has been computed for different energy bands, 
giving consistent results: the 0.3 — 1 keV and 3 — 10 keV en- 
ergy bands have been finally chosen to perform a one -to 
-one comparison to the results obtained by Margutti et al. 
(2010b) for flares detected in LGRBs. To this end, the LGRB 
flare time lags from Margutti et al. (2010b) have been re- 
calculated using the CCF analysis above (black dots in Fig. 
6): in Margutti et al. 2010b a pulse peak lag was instead cal- 
culated (grey dots in Fig. 6). The pulse peak lag is defined 
as — *peak - tpeak where *peak and tp 0ak are the peak 

times of the best fitting profiles in the energy bands I and 
II, respectively. As such, is sensitive to the assumed 

pulse fitting model: while the dependency is limited in cases 
of bright events, the limited statistics of the SGRB X-ray 
light-curves would cause the pulse peak lag computation to 
be inaccurate. For this reason we refer to the CCF time lag 
for both short and long GRB data, in the 7-ray and X-ray 
regimes. The light-curve time binning can potentially affect 
the derived for each SGRB the lag has been computed 
on light-curve pairs with 10 different time binnings spanning 
the range 0.2 — 20 s. The optimal time binning is defined as 
the lowest time scale giving origin to a CCF peak > 0.4 and 
is listed in Table 1. Larger binning time scales have been 
checked to produce consistent lag results. The window of 
time of investigation (ti and t/ of Table 1) has been de- 
termined selecting the time interval containing positive, at 
least 1 — a significant fluctuations around the smooth X-ray 
continuum (see Margutti et al. 2011 for details). 

For each SGRB, U and tf have been varied of ~ 20% 
both towards larger and lower values: consistent time lag 
values have been found. The sensitivity of the lag measure- 
ment to the smooth X-ray light-curve decay underlying the 
candidate flares has been investigated calculating the lag 
rf4 sub on light-curve pairs where the contribution of the 
smooth afterglow component has been properly subtracted 
and uncertainties propagated following the prescriptions by 
Margutti et al. (2011). For each SGRB this procedure has 
led to consistent r^ ub — values (r^ sub systematically 
has larger uncertainties due to the lower signal-to-noise of 
the subtracted light-curves). For this reason we refer to 
hereafter. Finally, we have tested and verified the robustness 
of our choice of energy bands to be compared, re- binning 
times and window of time of investigation, by performing a 
number of simulations where artificial lags have been first 
introduced into the light curves and then successfully recov- 
ered. Results axe reported in Table 1 . Figure 6 combines the 
T iag axid ^pk luminosity information in the lag-luminosity 
plane and clearly shows that SGRB flare candidates fall off 
the prediction based on LGRB flares and prompt pulses. 
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Figure 6. Lag-Luminosity plot. Red circles: OCF lag for candidate Hares of SGRBs: open symbols refer to SGRBs with detected EE; 
black (grey) triangles: CCF (pulse peak) lag for the sample of 9 flares of long GRBs of Margutti et al., 2010b, their Fig. 15. Black stars: 
prompt 7-ray data from the gold and silver sample of Ukwatta et al., 2010. Red squares: 3 a upper limits to the prompt lag of SGRBs 
for which it is possible to estimate the peak luminosity: open symbols refer to SGRBs with detected EE. The isotropic peak luminosity 
is computed in the 1 — 10 4 keV and 0.3 — lOkeV bands for prompt data and X-ray flares, respectively; the lag corresponds to the time 
difference between light-curve structures in the 50-100 keV and 100-200 keV channels (prompt data) and 0.3-1 keV and 3-10 keV channels 
(X-ray flares). Blue dashed line: best fitting law for the LGRB p\vmpt emission data. 


liable 1. X-ray and 7-ray properties for the sample of SGRBs analysed in this work. From left to right: GRB name: a (*) indicates an 
early-type host galaxy morphology (Fong et al., 2010; 2010b), while GRBs with detected extended emission are in boldface; redshift, 
duration and average spectral photon index of the prompt 15-150 keV emission as determined from GCNs; extended emission (EE) 
duration and EE to IPC (Initial Pulse Complex) intensity ratio (Hint) from Norris et al. 2010a: upper limits on i?i nt are listed when 
no EE has been detected (EE with zero duration); U and tf define the temporal window for the X-ray lag calculation while the At re b 
column reports the time scale used to re- bin the X-ray light-curve pairs; CCF time lag computed between 0.3 — IkeV and 3 — lOkeV; 
0.3 — lOkeV isotropic equivalent peak luminosity in the time interval U — tf as determined from the Norris 2005 profile fit (see Table 
Al); short lived (SL) or long lived (LL) X-ray afterglow according to the classification by Sakamoto & Gehrelg, 2009. 
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SL 

051227 


8.00 

1.31 ±0.22 

119.1 

0.0540 

101.2 

177.2 

4.00 

-4.5 ±7.2 

- 

LL 

060313 

- 

0.74 

0.71 ±0.07 

■ 0 

0.0005 

110.6 

250.00 

8.00 

30.5 ± 25.4 


LL 

070724A 

0.457 

0.40 

1.81 ±0.33 

0 

0.0074 

73.0 

126.0 

4.00 

3.6 ±4.6 

5.68 ± 2.32 

LL 

071227 

0.383 

1.80 

0.99 ± 0.22 

106.6 

0.0356 

126.18 

201.0 

3.00 

1.6 ±1.6 

1.13 ±0.44 

LL 

090607 


2.30 

1.25 ±0.30 

0 

0.0016 

76.1 

173.3 

10.0 

3.6 ± 10.4 


SL 

100117A* 

0.920 

0.30 

0.88 ±0.22 

0 

0.0030 

86.1 

238.5 

5.00 

3.3 ±1.6 

14.09 ± 5.60 

SL 


3 RESULTS 

The data analysis of the previous sections leads to the fol- 
lowing results: 

• SGRB flare candidates appear both in early-type and 
late- type host galaxy environments, irrespective of the short- 
lived (SL) or long-lived (LL) nature of the X-ray afterglow. 

• Both SGRBs with a bright extended emission (EE) and 
SGRBs which lack this component show cases of statistically 


significant fluctuations superimposed over smoothly decay- 
ing X-ray light-curves (Table 1). 

• Flares in LGRBs are known to show a spectral hard- 
ening during the rise time and a softening during the 
decay time, reminiscent of the prompt emission (e.g. 
Margutti et al. 2010b and references therein): as a result 
the hardness ratio (HR) evolution mimics the flare profile 
(see e.g. Goad et al. 2007, their Fig. 9) while the spectral 
photon index evolution anti- correlates with the flare flux. In 
spite of the lower statistics and limited A F/F of SGRB flare 
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Figure 7. Panel, (a): Relative variability flux A F/F vs. relative variability time-scale A t/t = w/t p k for the sample of X-ray flare 
candidates in SGRBs (filled circles), compared to early and late time LGRB X-ray flares (blue open diamonds and light-blue stars, 
respectively) from Chincarini et al., 2010 and Bernardini et al., 2011. The late time re-brightening detected in GRB 050724 is also shown 
for completeness with an orange open triangle (Bernardini et al., 2011). A small black dot marks data coming from SGRBs with detected 
extended emission. Solid, dashed and dot-dashed lines mark the kinematically allowed regions in different scenarios according to Ioka et 
al., 2005, their equations (7) and (A2). The A t/t and A F/F distributions are portrayed in panels (b) and (c) adopting the same colour 
coding. The green tick marks in panel (c) show the flux contrast for the sub-sample of LGRB flares of Fig. 6. 


candidates (Fig. 7), we find in the case with the best statis- 
tics a hint for correlation between the HR evolution and 
the temporal profile of GRB 101117A, with a photon index 
evolution that anti-correlates with the flux of the flare can- 
didates (Fig. 1, middle and lower panels). In the other cases 
(Fig. Bl, B2, B3, B4) the limited statistics prevents us from 
drawing firm conclusions. 

• The SGRB flare candidates width evolution is roughly 
linear in time and consistent with theu;/(l-}-z) vs. £ p k/(l+z) 
relation established by LGRB flares over ~ 4 decades in 
time (Fig. 2). The best-fitting law reads: tu/(l + z) = 
10~ 1 * o±o ‘ D (t pk /(l + 2 j)) 1 * 2±0 - 2 , It is remarkable that data 
coming from LGRB flares both at early and very late time 
(beyond i p k/(l + z) ~ 10 5 s, Fig. 2, inset) as well as tempo- 
ral fluctuations in completely different systems like SGRB 
afterglows are consistent at zero-order with the same, ap- 
proximately linear, law. We refer to Bernardini et al. 2011 
for a discussion of the possible bias affecting the w vs. t p k 
relation. 

* SGRB flare candidates are ~ 100 times dimmer than 


LGRB flares at the same rest frame time (Fig. 3). Selecting 
the sub-sample of LGRB flares detected in the same rest 
frame time interval 60s < t p k/(l + z) < 250s of SGRB flare 
candidates, we obtain a median (Lpkjp) ~ 10 49,8 ergs -1 to 
be compared to {L p k BB ) ~ 10 47,8 ergs -1 of the SGRB sam- 
ple showed in Fig. 4. As a result, SGRB flare candidates 
fall off of a factor ~ 100 the peak luminosity vs. time re- 
lation established by LGRB flares which reads: Lp^p B = 

l 0 5 4 - 8 ± 0 . 4( ipk /( 1 + ^- 2 . 7 ± 0.5 

• Short and long GRBs show a comparable 1 — 10 4 keV 
(rest frame) isotropic peak luminosity during their prompt 
emission, with a median (L p k,p) ~ 10 52,3 ergs -1 (Fig. 4, in- 
set). On the contrary, the peak luminosity of flares of both 
categories evaluated at the same rest frame time 60s < 
t P k/(l + z) < 200 s differ of a factor ~ 100 as noted above. 
While for LGRBs the typical prompt (1 — 10 4 keV) to flare 
(0.3-10 keV) peak luminosity ratio (Lp®p B /L B ^p B ) ~ 300, 
for SGRBs the same quantity reads: (L p k,p /£ p £f B ) ~ 
3 x 10 4 . Flare candidates in SGRBs are therefore less lu- 
minous than expected using the prompt-to-flare luminos- 
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ity scaling observed in LGRBs at the same flare rest frame 
time 5 6 * 

• SGRB flare candidates are ~ 100 times less energetic 

than LGRB flares (Fig. 5), with a median 0.3-10 keV en- 
ergy ~ 10 48 9 erg ({E\£% B ) ~ lO 50 9 erg). Since the 

width of LGRB flares and SGRB flare candidates are com- 
parable when evaluated at the same t p k/(l 4- z ), this result 
is a natural consequence of the (Lp£jP)/(Lpk f F B ) ~ 100 
reported above. On average, in the time interval 60s < 
t p k/(l + z) < 250s a flare in a LGRB emits ~ 0.6% of the 
1 - 10 4 keV prompt i?££p B ; the observed SGRB flare can- 
didates isotropic energy is ~ 0.2% the 1 — 10 4 keV prompt 
pSGRB 6 

‘E'iso.P 

• Flares and prompt pulses in LGRBs define a spectral 
lag-peak luminosity relation (Fig. 6): this finding is highly 
suggestive of a common origin (Margutti et al. 2010b) 7 . 
On the contrary, SGRB prompt pulses are known to ex- 
hibit much shorter lags than expected if they were to fol- 
low the LGRB prompt pulses lag-luminosity relation (e.g. 
Gehrels et al. 2006). Figure 6 extends this behaviour to their 
flare candidates: like SGRB prompt pulses, flare candidates 
in SGRBs fall off the lag-luminosity relation defined by 
LGRBs. 

• Flares candidates of SGRBs in the A F/F vs. rela- 
tive variability time-scale (A t/t = w/t p k) plane are com- 
patible with variability arising from density fluctuations of 
many regions viewed off-axis: on the contrary, neither the 
refreshed-shock nor the patchy-shell scenario is able to ac- 
count for the observed properties of the entire sample (see 
Ioka et al., 2005 for details on the definition of the various 
scenarios) 8 . In particular: a K-S test comparing the A t/t 
distributions of LGRB flares and SGRB flare candidates re- 
veals that they belong to the same parent population at 
~ 10% level of probability; the probability reaches the 88% 
level if LGRB flares are selected in the SGRB peak time 
range (60s < tpk/(l +z) < 250s). This is consistent with 
the common jj/(1 + z) vs. t p k/(l 4- z) relation followed by 
short and long GRBs discussed above. On the contrary, no 


5 The prompt emission peak luminosity is likely to be biassed 
towards the bright end of the real L p k,p distribution, since a min- 
imum signal- to- noise is needed to constrain the spectral parame- 
ters and calculate L p k, p in the 1-10 4 keV range: this requirement 
is more severe in the case of SGRBs whose observed emission is 
usually le^s bright than LGRBs. However, Fig. 8 shows that the 
main conclusion of this paragraph remains true even after relax- 
ing the requirement above. 

6 Again, this ratio is likely to be a lower limit to real value due 
to the bias affecting the the sample of SGRBs with prompt E- lso 
measure discussed in the previous paragraph. 

7 Note that flare and prompt lags are calculated in different en- 
ergy bands and are not directly comparable. The fact that flares 
define a lag-luminosity relation with slope very similar to the 
prompt data is however suggestive of a strict connection between 
flares and prompt pulses. See Margutti et al. 2010b for details. 

8 The smoking gun against a refreshed shock scenario would be 
the detection of a spectral change contemporaneous to the flare 
candidates: while the statistics of the XRT light-curves of SGRBs 
is limited in the case with best statistics the HR evolution is 
correlated to the candidate flare profiles and the continuum after 
the flaring emission is softer than the emission detected during 
the period of temporal variability (Fig. 1). These findings favour 
an alternative explanation. 


Table 2. Summary of the properties of SGRB flare candidates 
compared to LGRB X-ray flares observed at the same t p k/(14^). 


Property 

= 

* 

Width 

Relative variability 
Relative variability flux 

w(tpk) 
A t/t 

AF/F 

Peak luminosity 


t SGRB 

~ 0.01 

pk,F 

Isotropic energy 

Lag-luminosity 

Flare to prompt energy ratio 

■Etso.F 

E iso,P 

piSGRB 
Lpk (^lag) 

Flare to prompt luminosity ratio 
Flare to prompt pulse width ratio 


Lpk,F/Lpk,P 

Wp/wp 


X-ray flare candidate in a SGRB shows a relative variability 
flux A F/F > 2 in strong contrast with the LGRB A F/F 
distribution at comparable i p k/(l + z) (Fig. 7). Such promi- 
nent flares would be easier to detect, so that it is unlikely 
that an observational bias could explain the present lack 
of detection. A K-S test comparing the two distributions 
shows that the probability that LGRB and SGRB flare can- 
didates share the same A F/F parent population is as low as 
3.3 x 1CT 4 . This result partially inherits the uncertainty af- 
fecting the completeness of both distributions for very small 
A F/F values. Another source of uncertainty arises from the 
difficulty in evaluating the continuum underlying the po&- 
sible flare emission in SGRBs when data are particularly 
sparse. In spite of these limitations, after more than 6 years 
of Swift observations (and ~ 60 SGRB afterglows detected) 
there is still no SGRB showing a prominent (A F/F > 10) 
fast-varaibility A t/t 1 feature during its X-ray afterglow. 

The SGRB flare candidate A F/F is instead more similar to 
the relative variability flux of flares in LGRBs detected at 
late times (t p k > 1 ks, light-blue stars in Fig. 6, main panel; 
Bernardini et al. 2011): the two A F/F distributions share 
the same parent distribution at ~ 21% level of probability 
(K-S test). 

The results above demonstrate the complexity charac- 
terising the SGRB flare candidates phenomenology: Table 
2 reports a summary of their properties when compared 
to LGRB flares. The main result is that the population of 
SGRB X-ray flare candidates only partially share the set 
of observational properties of LGRB X-ray flares detected 
at the same rest frame time: are there real X-ray flares in 
SGRBs? A detailed discussion is provided below. 


4 DISCUSSION 

Observations show that like LGRBs, at least some SGRB X- 
ray afterglows deviate from a smooth power-law decay and 
show variability. In the following we discuss the properties 
of SGRB flare candidates providing a one-to-one comparison 
with LGRB X-ray flares. The aim is to better understand the 
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origin of the short burst afterglow variability and uncover 
potential links with the prompt phase. 

With A t/t > 1, the late-time (£ p k ~ 5 x 10 4 s) 
re-brightening of GRB 050724 strongly differs from the 
properties of the entire sample of flare candidates de- 
tected in SGRBs thus questioning its classification as 
flare-like episode (Panaitescu et al. 2006: see however 
Grupe et al. 2006; Campana et al. 2006; Malesani et al. 
2007). Malesani et al. (2007) report the detection of an 
optical and radio re^brightning associated to the X-ray 
bump which is unusual if compared to the standard prop- 
erties of X-ray dares, while being more common to late- 
time re-brightneings observed in LGRBs as well (see e.g. 
GRB 081028, Margutti et al. 2010a). In addition, no hard- 
to-soft evolution can be detected in the X-ray data (Evans 
2010), which is instead typical of flares (Goad et al. 2007) 
and prompt pulses in LGRBs (Hakkila et al, 2011). In the 
following we focus our attention on SGRBs fast variabil- 
ity (A t/t < 1) referring the reader to Bernardini et al. 
(2011) for a complete discussion of the late-time behaviour 
of GRB 050724. 


4.1 SGRB local and global environment 

The standard model (see Nakar et al. 2007 for a recent re- 
view) explains the X-ray afterglow of long and short GRBs 
as synchrotron radiation arising from the deceleration of 
a relativistic blast wave into the external medium. If the 
shock front is homogeneous and expands into a smooth am- 
bient density, a smooth afterglow light-curve is expected. 
In this context, variability in the X-ray afterglow could be 
caused by re-freshed shocks (i.e. shocks caused by slow shells 
catching up with the leading, decelerated shell at late times, 
Kumar k Piran 2000a, Granot al. 2003): Fig. 7 shows that 
half of the SGRB flare candidates sample do not comply with 
this scenario 9 . Furthermore, the spectral variability shown 
in Fig. 1‘, lower panel, makes it difficult to interpret the flare 
candidates in the re-freshed shocks scenario. 

A first alternative is to relax the assumption on the 
homogeneity of the shock front (Kumar k Piran 2000b): an 
intrinsic angular structure of the emitting surface is able to 
produce variability with a characteristic time-scale At ^ t 
if the angular fluctuation is persistent (patchy shell model, 
Nakar et al. 2004) 10 : no SGRB flare candidate is consistent 
with this expectation (Fig. 7, main panel). 

A second alternative invokes the presence of ambient 
density fluctuations either caused by turbulence in the ISM 
or by variable winds from the progenitor. From Table 1 
it is however apparent that temporal X-ray variability has 
been detected for SGRBs residing both in early-type and 
late-type host galaxies which likely have very different ISM 
properties. In particular GRB 050724 and GRB 100117A are 
the unique two SGRBs with secure early-type host galaxy 
association (Barthelmy et al. 2005b; Fong et al. 2010). This 


9 See however Granot al. (2003); Vlasis et al. (2011): sharp op- 
tical and radio flares could be produced by collision of ultrar 
relativistic shells. 

10 Details on variability arising from a time- varying anisotropic 
emitting surface can be found in Ioka et al., (2005). 


suggests that the ISM turbulence is unlikely to provide a fear- 
sible physical mechanism for the detected variability. Note 
however that the limited size of our sample prevents us 
from quantitatively discussing the correlation between the 
appearance of flare candidates and host environment. 

Different progenitor models of SGRBs lead to distinct 
predictions on their local environment as well. In particu- 
lar, according to the standard compact binary merger sce- 
nario (NS-NS or NS-BH, Eichler et al. 1989; Narayan et al. 
1992), no variable wind is expected from the progenitor. 
An origin of SGRB flare candidates from density fluctua- 
tions of the circumburst environment is therefore consid- 
ered unlikely. Alternative scenarios leading to SGRB envi- 
ronments with circumburst densities comparable to LGRBs 
have been however explored by Nysewander et al. (2009) to 
explain the similar Fr/Fx ratio (where Fr and Fx stand 
for the afterglow flux density at 11 hrs post-trigger in the 
R and X-ray bands, respectively). A systematic difference 
between sub-galactic environments able or unable to pro- 
duce variability in the X-ray afterglow could in principle 
be revealed by different offsets with respect to the host 
galaxy centers: however, the present sample of SGRBs with 
flare candidates includes both SGRB with large offsets (e.g. 
6 = 14.80 ± 0.34 kpc for SGRB 071227, Fong et al. 2010b; 
D’Avanzo et al. 2009) and events with very small offsets 
(e.g. S = 0.47± 0.31 kpc, SGRB100117A, Fong et al. 2010). 
While the observed offset distribution is incomplete, these 
data suggest that the properties of the local environment of 
SGRBs are not the key parameters determining the presence 
of flare candidates in their X-ray light-curves. This conclu- 
sion is strengthened by the results of Sect. 4.5. 

Flares candidates and prompt pulses could alternatively 
share a common origin. In this case, flare candidates would 
bring no information on the external medium density. We re- 
fer to this interpretation as the ” internal origin” possibility 
(see Sect. 4.7). 

4.2 Flare candidates and extended emission 

3 bursts in our sample (boldface in Table 1) present evidence 
for extended emission (EE): a long-lasting (At ~ 10 2 s) soft 
X-ray tail that follows the short hard spike in the prompt 
phase. Norris et al. (2010a) analysed a sample of 51 Swift- 
BAT SGRBs looking for the presence of EE in their 7-ray 
data and concluded that ~ 3/4 of the BAT SGRBs are not 
accompanied bv an EE component. In particular, in this 
work it is shown that the upper limit on the EE to IPC 
(initial pulse complex) intensity ratio of the SGRBs of our 
sample (Le. with flare candidates) not showing EE is a factor 
> 10 below the standard values found for GRBs with EE 
(Norris et al. 2010a, their Fig. 1). This implies that in those 
cases the EE is either very weak or absent. The inverse is 
also true: some SGRBs with bright EE are accompanied by a 
smooth X-ray light-curve at t > 80 s (see e. g. GRB 080503, 
Norris et al. 2010a, their Table 1). 

While it is still unclear if the origin of the EE compo- 
nent is related to the prompt emission, the afterglow or a 
third, unknown physical process (Perley 2009; Norris et al. 
2010a; Norris et al. 2011), our analysis shows that the pres- 
ence of a bright EE does not imply the presence of flare 
candidates and viceversa. We stress that the higher aver- 
age brightness characterising the XRT light-curves of GRBs 
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with EE with respect to those without EE, would naturally 
bias the result in the opposite direction (i.e. it would favour 
temporal fluctuations to be detected in XRT light-curves 
with the best statistics, leading to a biassed flare-EE con- 
nection) thus strengthening our conclusion. The limited size 
of our sample (which is however the widest possible at the 
time of writing) does not allow us to quantitatively discuss 
the flare candidates vs. EE correlation (or lack thereof). 

Our analysis cannot however exclude that SGRB flare 
candidates constitute temporal fluctuations superimposed 
on (and physically linked to) the X-ray tail of the EE, 
as suggested by the epoch of flare candidates detection 
(t p k/(l + z) < 200 s). Furthermore, the limited range 
spanned by the flare-to-continuum flux ratio A F/F (with 
A F/F ~ 1, Fig. 7) is suggestive of a physical link between 
flare candidates and the underlying continuum (i.e. the EE). 
Again, flare candidates would be associated to both bright 
and faint (undetected in the 7-rays) EEs. 

4.3 Time scales 

The relative time-scale distribution A t/t of SGRB flare 
candidates (Fig. 7, upper panel) is compatible with being 
drawn from the same parent population of flares detected 
in LGRBs at similar t p k/(l + 2) at ~ 88% level of proba- 
bility: with a median value of (A t/t) = 0.25 and extending 
from 0.1 up to 0.5 it furthermore satisfies the limit A t/t ^ 1 
which is the classical argument against fluctuations arising 
from the external shock (e.g. Zhang et al. 2006). The A t/t 
distribution therefore does not support an external shock 
origin for flare candidates in SGRBs (see however Dermer 
2008), in agreement with the findings of Sec. 4.1. 

The evolution of the flare duration with time w(t p k) for 
t P k > T$c (Fig. 2) observationally distinguishes flares from 
prompt pulses in LGRBs (Margutti et al. 2010b): as time 
proceeds, LGRB flares becomes wider and wider following 
the best fitting law: w/(l + z) ~ (t p k/(l + ^)) 1 ‘ 2±0 ’ 2 . This 
quasi-linear regime establishes for £ p k/(H-z) > 290/(1 + 2), 
and extends 4 decades in time up to t P k / ( 1 + 2) > 10 5 s. No- 
tably, while the i^ G p B released during the LGRB prompt 
emission spans more than 3 decades in energy (Fig. 5), likely 
reflecting different properties of the LGRB central engines, 
the duration of subsequent episodes of activity seems to fol- 
low a universal evolution with limited dependence on the 
initial energv budget. Flares with different amplitudes A, 
flux ratios A F/F and fluences show similar w/(l + 2) when 
appearing at the same t p k/(l + 2) (Chincarini et al. 2010), 
suggesting that the main parameter driving the flare width 
evolution in a LGRB is the elapsed time from the explosion 
t P k/(l + 2) 11 . 

Figure 2 shows that flare candidates in SGRBs are con- 
sistent with the same quasi-linear temporal scaling: from the 
width measurement it is not possible to distinguish a flare 
in a LGRB from a flare candidate in a SGRB. The tem- 
poral properties of the prompt emission of long and short 
GRBs are clearly different in terms of duration and vari- 
ability (e.g. Nakar et al. 2002): however, ~ 30 s later, the 

11 As a second level of approximation, one should consider that 
more prominent flares (A 1) tend to be wider (Margutti et al. 
2010b). 


width of fluctuations superimposed over their X-ray light- 
curves seems to have lost any information on the nature of 
central engines able to produce 7-ray photons with such dif- 
ferent temporal properties. In both cases the w evolution is 
driven by the £ p k, irrespective of their different initial con- 
ditions (and initial variability time scales): while for flares 
wf GRB /(l + z) ~ typ GRB /(l + z) at similar t p k /(l +2), 
for prompt pulses u>p GRB < wp GRB (with Wp GRB ~ 0.05 
s to be compared to wp GRB ~ 1 s, observer frame values, 
Nakar et al. 2002; Nakar et al. 2002b 12 ). 

This finding strongly suggests that the origin of the 
quasi-linear w/( 1 + z) vs. t P k/(l + 2) scaling must be within 
what is in common for the long and short GRB models, 
irrespective of the progenitors, environment, lifetime, vari- 
ability time-scales and prompt energy release of their cen- 
tral engines. Both models require the prompt emission to 
originate from a relativistic outflow (see e.g. Piran 2004): 
if the longer and longer duration of flares is of dynamical 
origin and dominated by the expansion of the emitting re- 
gions, no memory of the properties of the central engine 
which launched the relativistic outflow would be preserved 
so that long and short GRBs would display flares with sim- 
ilar u;(£ p k). Alternatively, in the context of accretion mod- 
els the iu(i p k) relation originates from the viscous evolution 
of the hyper-accreting disk around the newly born black 
hole, common ingredient of the likely progenitors of the two 
classes (Perna et al. 2006). Our results would imply a strict 
analogy between the mechanisms that regulate the late-time 
evolution of accreting disks originating from collapsars and 
mergers of compact objects, irrespective of their different 
masses and life times (according to the standard scenario 
accreting disks related to SGRBs are likely to be less mas- 
sive then LGRB disks and short-lived, as suggested by the 
observed prompt duration. See Nakar et al. 2007 for details). 

Lazzati et al. (2011) recently suggested that instabili- 
ties arising from the propagation of the jet through the dis- 
rupting star could explain the presence of flares in LGRBs 
with A F/F < 10. Even assuming a continuous and feature- 
less release of energy by the central engine, the pressure of 
the surrounding stellar material would provide the physical 
origin for the intermittent flare behaviour, naturally explain- 
ing the universal (i.e. with limited sensitivity on the star 
properties and energy budget) quasi-linear w(t p k) relation 
(Lazzati et al. 2011, their Eq. (6)). However this model can- 
not account neither for presence of flares in SGRBs nor for 
their similarity to LGRB flares in the w/(l+z) vs. t p k/(l+2) 
plane: according to the majority of SGRB progenitor scenar- 
ios the engine is ” exposed 5 and there is no stellar material 
the jet could interact with. As a consequence, if the iu(t p k) 
relation in LGRBs originates from the interaction with the 
progenitor stellar material, it is difficult to explain why flare 
candidates in SGRBs, while originating from a completely 
different physical scenario, are however consistent with the 
same relation, as observed. Our results therefore imply that 
either the LGRB flare w(t p k) relation does not originate 
from the jet-stellar material interaction or that the progen- 


12 Note however that the first 1-2 s of bright LGRBs display 
variability time scales comparable to the SGRB prompt emission, 
Nakar et al. (2002). 
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itors of long and short GRBs are much more similar than 
previously thought (see e.g. Lazzati et al. 2010). 

Finally, our findings suggest that while the variability 
time scales measured during the prompt emission could still 
directly reflect the intrinsic variability of the central engine 
(see Piran 2004 and references therein), absolute measures 
(i.e. not re- normalised) of flare widths likely do not (if time- 
dilated by physical mechanisms which are only indirectly 
related to the central engine). On the contrary, the ratio of 
interesting time-scales of the same flare profile, being sub- 
ject to the same temporal stretching, could still retain an 
imprint cf the original mechanism at work: this would ex- 
plain why in LGRBs the flare rise time t T and decay time 
like the w , linearly evolve with t p k (Chincarini et al. 2010; 
Bernardini et al. 2011) preserving their ratio t r ~ 0.5£<j over 
4 decades in £ P k and leading to flares with asymmetry val- 
ues very similar to the prompt pulses (while being a factor 
^ 100 wider) 13 . 

4.4 Flux contrast 

While the w(f P k) and A t/t measurements do not allow us 
to distinguish a flare in a LGRB from a flare candidate 
in a SGRB 14 , the flux contrast distributions A F/F of the 
two populations are significantly different (Fig. 7), with 
SGRB flare candidates having systematically lower AF/F 
values. In contrast to LGRBs, none of the ~ 60 Swift X- 
ray afterglows associated to SGRBs shows cases of powerful 
(A F/F » 5), highly variable A t/t < 1 flares. A Chandra 
observation of 9 X-ray photons from SGRB050709 15 ~ 15 
days after the explosion led Fox et al. (2005) to conclude 
the presence of high-amplitude (A F/F « 10), fast variabil- 
ity (A t/t « 0.01) in its X-ray afterglow. Swift observations 
suggest that this kind of variability, if present, is extremely 
rare. SGRBs lack the presence of prominent fast- rise expo- 
nential decay features superimposed over their X-ray after- 
glow for : p k/(l + z) > 60 s. 

Lazzati et al. (2011) predicted that if flares in SGRBs 
originate from the intrinsic variability of their inner engine, 
their AF/F distribution should be analogous to the brighter 
population (A F/F ~ 100) of LGRB flares. With a maxi- 
mum A F/F « 2 the detected SGRB flare candidates are 
not consistent with these expectations and populate instead 
the low end of the LGRB distribution unless the SGRB con- 
tinuum flux F has been overestimated by a factor ~ 100 
which we consider unlikely. An interesting possibility is that 
the X-ray light-curves of SGRBs are dominated by an emis- 
sion component which is not present in the LGRB after- 
glows (such as the EE): this would lead to systematically 
lower A F/F for SGRBs when compared to LGRB flares. 

13 Note that given the limited statistics of the SGRB X-ray af- 
terglov’s, nothing can be said about the asymmetry of SGRB flare 
candidates. 

14 The limited number of flare candidates plays a role in this 
statement. We cannot exclude that a significant improvement of 
the SGRE and LGRB statistics could lead to the detection of 
secondary effects. 

15 GRB 050709 is a short burst detected by HETE for which 
Swift did not do the follow up. For this reason it is not included 
in the present sample. The episode of prominent variability is 
possibly related to a statistical fluctuation. 



Figure 9. Relative variability flux A F/F evolution as a func- 
tion of ipk/Tgo: black triangles and green stars indicate LGRB 
flares with t p k < 1 ks (Chincarini et al., 2010) and t p y > 1 ks 
(Bernardini et al., 2011), respectively. Red open and filled circles: 
flare candidates in SGRBs with and without extended emission, 
respectively. Horizontal dashed, dot and dot-dashed lines mark 
the A F/F median value for the three samples. 

The clustering of the A F/F ratio around 1 would suggest 
a link between flare candidates and underlying continuum. 
However, the strong correlation found by Nysewander et al. 
(2009) between the X-ray afterglow flux and prompt fluence 
of both long and short GRBs reveals a high degree of similar- 
ity between their X-ray afterglows at least at late times (i.e. 
11 hrs, observer frame). A detailed comparison of short and 
long GRB X-ray afterglows at early times is in preparation. 

The SGRB A F/F observed at £ p k ~ 100 s is instead 
typical of LGRB flares detected at much later times: tpk > 
Iks (Bernardini et al. 2011). From Fig. 9 it is apparent 
that SGRB flare candidates are consistent with the A F/F 
distribution of LGRB flares detected at the same £ p k/Tgo 
epoch (where Tgo is the duration of the prompt 15-150 keV 
emission): differences instead arise if we compare the prop- 
erties of the two classes at the same £ pk /(l 4- z). While 
the T90 is possibly a crude estimation of the intrinsic time 
scale of evolution of the central engine, this result opens 
the possibility that prominent fluctuations are not currently 
detected in SGRB afterglows because of the intrinsically 
late-time Swift re-pointings: £ r ep ^ 100 s corresponds to 
trep/Tgo ~ 100 - 1000 (SGRBs) and to t rGp /T 90 ~ 1 - 10 
(LGRBs). From another perspective this finding directly 
connects the flux properties of SGRB flare candidates to 
the evolution time scale of their central engines. The prompt 
Tgo qualifies as a good proxy for the intrinsic time scale that 
drives the subsequent flux evolution of flares and flare candi- 
dates in long and short GRBs, respectively. This conclusion 
is strengthened by the results of Sec. 4.5. 

4.5 Energetics 

4 5.1 Flare L p k evolution with time 

While SGRB prompt pulses compete with those of LGRBs 
in terms of peak luminosity, the same is not true for their 
late- time variability: SGRB flare candidates are a factor 
~ 10 2 dimmer that expected (Fig. 3). Sec. 4.4 and Fig. 

5 suggest that long and short GRBs might be astrophysi- 



12 R. Margutti et al. 



Figure 8. Re-normalized peak luminosity vs. re-normalized peak time for the sample of LGRB flares from Chincarini et al. (2010) and 
Bernardini et al, (2011) (triangles) with ipk/Tbo 5 10 3 and SGRB flare candidates (circles). Dark and light colours distinguish events 
with and without redshift measurement to allow a direct comparison with Fig. 3. Open symbols refer to SGRBs with detected extended 
emission. Black dot dashed line: best fitting power-law model for the LGRB sample: L p k,F/(Lp) = lO -0 ' 72 ^ 0 ' 03 ^^/^^) -0,9 ^ 0 ’ 1 and 
extrinsic scatter a — 0.62 ±0.01. The shaded area marks the ±1 a region around the best fit. From top to bottom, the blue dotted lines 
mark the ^iso,F/^iso,P = 1,0.1,0.01,0.001 regions of the plane, where Fliso.F and F?i so ,p have been calculated in the 0.3-10 keV and 
15-150 keV energy bands, respectively. The flare width vs. peak time relation from Chincarini et al. (2010) has been used. 


cal systems evolving on completely different time-scales and 
with different energy budgets but otherwise based on a sim- 
ilar physical mechanism of emission. In that case, we would 
expect the evolution of the flare luminosity to exhibit much 
better agreement between long and short GRBs once the in- 
trinsic time-scale and central engine energy scaling are prop- 
erly accounted for. Figure 8 shows that this is indeed the 
case: long and short GRBs are consistent with a common (al- 
beit highly scattered) L p k,F/(Lp) ~ (fpk/Tbo) -0,9 * 0,1 scal- 
ing (where (Lp) and T90 are the average prompt luminosity 
and duration of each GRB in the 15-150 keV energy range, 
respectively). Figure 8 therefore establishes a direct connec- 
tion between the properties of SGRB flare candidates and 
LGRB flares, providing further support to a common, likely 
internal origin (see Section 4.7). 

In the context of a one to one comparison of short vs. 
long GRB X-ray variability, the re-normalised peak lumi- 
nosity and re-normalised peak time are the most relevant 
quantities: the observed T90 and (Lp), with their sensitivity 
to the instrument threshold and energy band passes, likely 
provide crude (but nevertheless the best) approximations 
to the exact values to be used (while being partially re- 
sponsible for the large scatter of the relation). In particu- 
lar the loss of total fluence in the prompt 7-ray due to the 
limited Swift - BAT band-pass likely affects short more than 
long GRBs. As a result the SGRB L p k,F/(Lp) values may be 
overestimated when compared to LGRB values 16 . However, 

16 This effects would mainly depend on the different prompt 


Fig. 8 shows that even a factor of several of overestimation 
(Nysewander et al. 2009) would not strongly affect our con- 
clusions. 


4 - 5.2 Implications of the L p k,F/(Lp) vs. tpk/Tgo relation 

SGRB flare candidates are consistent with the highly scat- 
tered L p k,F/(Lp) vs. tpk/Tgo relation established by LGRB 
flares: since the origin of LGRB flares is likely connected 
to their prompt pulses (Krimm et al. 2007; Margutti et al. 
2010b), we speculate that a similar physical mechanism (ex- 
cept for its energy budget and life-time) powers long and 
short GRBs: observation ally speaking, the main distinction 
between flares and prompt emission in LGRBs is the evo- 
lution of the former with time for t p k > Tqq. It is therefore 
natural to expect a similar scaling of LGRB flares and SGRB 
flare candidates in terms of tpk/Tbo if they share a common 
origin. 

A comparison of the prompt 7-ray emission of SGRBs 
to the initial 2 s of LGRBs reveals a high degree of sim- 
ilarity in the pulse duration distributions and correlations 
in the temporal structure of the two classes (Nakar et al. 


emission spectral peak energy £p e ak of long and short GRBs: 
however, Ghirlanda et al. (2009) showed that SGRBs have a 
harder low-energy spectral component but only slightly higher 
when compared to LGRBs. 
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2002; McBreen et al. 2001) 17 ; an analogous study was per- 
formed by Ghirlanda et al. (2009) in the spectral domain: 
based on the spectral analysis of the prompt emission of 79 
short and 79 long GEBs detected by BATSE the authors 
showed that no difference is found comparing the spectral 
properties of SGRBs to the first 1-2 s emission of LGRBs. 
Temporal and spectral studies therefore point to a common 
mechanism operating in the first few seconds of any event. 
The present work extends this similarity to their late-time 
emission. 

The quasi-linear ui(t p k) evolution shared by long and 
short GRBs allows us to draw reference lines of equal flare- 
to-prompt energy ratios as a function of t p k/79o in the 
L p ^ t p/(Lp) plane (blue dotted lines of Fig. 8). Flares in the 
0.3-10 keV band pass are found to emit between (0.1—100)% 
of the prompt 15-150 keV Eiso.p, with the majority of them 
lying between the 1% and 10% levels. Little evolution of 
the flare-to-prompt energy ratio E' 1SOi f / E i so ,p in terms of 
tpk/Tvo can be inferred from the plot (for tpk/lho < 300). 
In particular, SGRB flare candidates and LGRB flares show 
comparable flare-to-prompt energy ratios, as reported in Ta- 
ble 2. This finding provides further support to a physical link 
between LGRB flares and SGRB flare candidates. 

The GRB central engines seem to release comparable 
fractions of prompt emission energy at late times, irrespec- 
tive of the long or short GRB nature. From Fig. 8: 


L p k,F 

(L P ) 



( 1 ) 


with a = —0.1 ± 0.1. The best-fitting rest-frame w(t p k) re- 
lation reads: 


(iSiHGfc) 


0»+i) 


( 2 ) 


with fi = 0.2 ± 0.2 (Fig. 2). Equations (1) and (2) express 
common properties of LGRB flares and SGRB flare candi- 
dates. The normalisation parameters iVi and N% possibly 
hide the dependence of L P k,F/(Lp) and w on other parame- 
ters. This hidden dependence might be partially responsible 
for the large scatter of relation (1); ot and ft parametrise 
the non-linear dependence of Lpk f F/{Lp) and w/{ 1 -I- z) on 
tpk/Tgo and t p k/(l-|-z), respectively. Combining Eq. (1) and 
Eq. (2) it is possible to show that: 


ffiso,F 

Eiso,P 



ot 

(3) 


where N ~ 0.9 2 iViAr 2 and .Eiso.F ~ 0.9L p kU>/(l + z) (valid 
for a Norris et al. 2005 flare profile where w is calculated be- 
tween 1/e intensity points and with £ r = 0.5id as found by 
Chincarini et al. 2010 for LGRB flares). Equation (3) shows 
that Ei SOi F/Ei SOf p oc [t P k/(l + a weak dependence of 

the flare-to-prompt energy ratio on t P k/(l 4- z) is expected 
if P — a rsj 0. From the best-fitting relations, we find that 
both parameters are consistent with 0 at lcr: in particu- 
lar 0 — a = 0.3 ± 0.2 (consistent with 0 at 1.5 a level). At 
similar i P k/(l + z) a residual dependence on the rest-frame 
prompt duration Tqo/(1 4- z) is expected to arise from the 


17 The. temporal evolution of pulses as a function of frequency 
(i.e. the spectral lag) shows however dissimilarities, as discussed 
in Section 4.6. 


third term of Eq. (3): this dependence, if present, would 
be able to distinguish the population of LGRB flares from 
SGRB flare candidates detected at the same rest-frame peak 
time in terms of Ei SOl p / E] so ,p ■ However, a = —0.1 ±0.1 (and 
the relation is highly dispersed). We therefore conclude that 
the quasi-linear w{t p k) and L p k,F/(Lp) vs. £ p k /Tqo relations 
translate into i?iso,F/-Eiso,p ratios which, at first order ap- 
proximation, show limited dependence on both the proper- 
ties of the central engine (i.e. duration of the prompt emis- 
sion) and elapsed time from the explosion. 

From Fig. 3, L p k,F B /^pk^F B ~ 10“ 2 at the same t P k/(l4- 
z). A factor ~ 100 is roughly the ratio of the isotropic energy 
emitted by long and short GRBs during their -prompt 7-ray 
emission (Fig. 5, inset) 18 . LGRB flares and SGRB flare can- 
didates are therefore expected to show a similar behaviour 
in the L p k,F / E' lS0f p vs. t p k/(l + z) plane (in strict analogy 
with the X-ray afterglow scaling found by Nyse wander et al. 
2009, their Fig. 6). Equation (1) can be easily re-arranged 
into: 


w =N f ^ y (a+1 Y r 9 Q y 

tfiso.p Vi+zy ’ 


(4) 


Again, the limited departure of L p k,F/{±p) from a linear 
relation in t p k/Tgo (a = —0.1 ±0.1) causes LGRB flares and 
SGRB flare candidates to share the same scaling (at least 
at first order approximation) . 

All the above indications point to a SGRB flare candi- 
dates internal origin (Section 4.7) and establish a connection 
between long and short GRB X-ray variability. 


4.6 The lag- luminosity relation 

Negligible spectral lag above ~ 25 keV is the fundamental 
attribute that makes the prompt 7-ray emission of short 
bursts different from LGRBs, in addition to their nar- 
rower pulses, shorter duration and slightly harder emission 
(Norris & Bonnell 2006 and references therein). The spec- 
tral lag is the time difference between the arrival of high- 
energy and low-energy photons: in our analysis, a positive 
value indicates that high energy photons lead the low en- 
ergy emission. During the prompt phase of LGRBs the spec- 
tral lag Tj^ is anti- correlated with the peak luminosity as 
shown by Fig. 6 (Norris et al. 2000; Ukwatta et al. 2010); 
in contrast, short bursts have small rj^ (Norris & Bonnell 
2006) and occupy a separate area of the L p k vs. t^_/(1 + z) 
parameter space (Fig. 6 and Gehrels et al. 2006). Recently, 
Margutti et al. (2010b) have demonstrated that, in strict 
analogy to their prompt pulses, LGRB X-ray flares define a 
lag-luminosity anti-correlation, where the lag is computed in 
the X-ray band (black dots of Fig. 6). With the present work 
we complete the observational picture above, showing that 
flare candidates in SGRBs fall off the lag-luminosity relation 
defined by LGRBs: this furthermore supports a robust con- 
nection between prompt pulses and flare candidates in short 


18 Note that the same ~ 100 factor is found as the ratio of the 
isotropic energy emitted by LGRB flares and SGRB flare candi- 
dates: this is a direct consequence of their flare peak luminosity 
ratio ~ 10^ 2 , coupled to the very similar w(t p k) 

evolution. 
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bursts. At the same time this result points to some differ- 
ences between LGRB flaxes and SGRB flare candidates. 

While SGRB prompt pulses are significantly narrower 
than LGRB pulses (Nakar et al. 2002; Nakar et al. 2002b), 
flares show instead comparable width in both classes {Fig. 
2). Hakkila et al. (2008) showed the existence of a lag- width 
correlation for prompt pulses: the wider the prompt pulse, 
the longer the lag. This behaviour has been recently ex- 
tended to LGRB flares by Margutti et al. (2010b). The sim- 
ilar width of LGRB flares and SGRB flare candidates im- 
plies that the lag- width relation cannot be invoked to explain 
the lag of SGRB flare candidates which are systematically 
shorter than expected from the lag- luminosity relation of 
LGRB flares. 

The physical cause of lags in the GRBs prompt 
emission is not yet understood: variations in the line-of- 
sight (Salmonson 2000); variations of the off-axis angle 
(Ioka & Nakamura 2001) and rapid radiation cooling effects 
(Schaefer 2004) are a few of the proposed models. In partic- 
ular it is at the moment unclear if lags in short bursts are 
small and non-measurable or intrinsically zero. According to 
the first possibility short and long GRBs would be powered 
by a similar, progenitor-independent physical mechanism, 
with SGRBs being faster evolving versions of LGRBs. The 
latter would instead point to some intrinsic differences. At 
present it is not possible to observationally distinguish be- 
tween the two scenarios 19 . 

In the case of flares, the situation is complicated by 
the limited and fixed energy band-passes used for the lag 
calculation (0.3-1 keV vs. 3-10 keV). The fundamental ori- 
gin of the observed lag is the spectral evolution of a pulse 
(or flare) profile to lower energies (Kocevski &; Liang 2003; 
Margutti et al. 2010b). The spectral peak energy E pea k de- 
crease in ume plays a major role in determining the observed 
spectral evolution and lag value (see Margutti et al. 2010b 
for details). As a consequence, if the observed J5 pea k does 
not cross the instrumental band pass during the emission, a 
limited spectral evolution will be detected and a short time 
lag determined. The shorter (when compared to what ex- 
pected from the LGRB flare lag-luminosity relation) time 
lag of SGRB flare candidates might be partially a conse- 
quence of this observational effect (while possibly being in- 
trinsically larger) 20 : the results from Section 4.4 and 4.5 
suggest that long and short bursts basically differ in the in- 
trinsic time scale of central engine evolution (with SGRBs 
evolving faster). Since for LGRBs ISpSJV^p^k^ iS 0-01 
(with £p^k,F ~ 1 — 3 keV, observed value, Margutti et al. 
2010b), the faster evolution of short bursts likely implies 
-® P Sk,F/^pS.k^ < 0.01 for flares detected at the same ob- 
served This result translates into: £p^k ( F < 1 keV con- 


19 FYom the prompt lag-v/idth relation of Hakkila et al. (2008): 
r 7 , LGRB LGRB ™ 0.01 — G.l, Using u;| GRB ~ 0.05 s as typical 
value from Nakar et al. (2002b), we have t 7 iSGRB ^ q . 5 — 5 ms 
assuming that the r' 1 /w p ratio is universal. We typically resolve 
lags in SGRBs with a sensitivity of a few ms. This implies that 
we would be barely able to measure lags in the SGRB prompt 
pulses even if SGRBs were to follow the LGRB t"* jw p relation. 

20 Note that the limited brightness of the flare candidates com- 
pared to the underlying X-ray continuum does not allow us to 
perform a one-to-one comparison with the spectral properties of 
LGRB flares to quantitatively check this possibility. 


sidering that ~ 33 OT< ^ CT of magnitude esti- 

mation (Ghirlanda et al. 2009). According to this scenario, 
^pSk,F ^ below the XRT band for the majority of the emis- 
sion, possibly leading to a lag underestimation. The presence 
of this observational bias makes the interpretation of the en- 
tire lag-luminosity relation far from being straightforward. 
We stress that the dependence of the lag-luminosity on the 
choice of the fixed energy bands (both in the rest frame and 
in the observer frame) should be removed before addressing 
the physical interpretation of the anti-correlation. However, 
this topic goes beyond the scope of this paper and will be 
addressed in a future work. 

With this caveat in mind we note that if the energy 
JEJiso.F released by flares at different i p k is similar (as in- 
dicated by Section 4.5), then, considering that L p k,F ~ 
Ei SO ,F/w with the lag being positively correlated to the w 
(Margutti et al. 2010b), we would expect L p k,F to be anti- 
correlated with 7i* g as observed for LGRB flares of Fig. 
6: L p k,F = A^^riag) -1 . In particular the normalisation 
iViag oc ISiso.F, which implies g RB ~ lOOiV^g 118 (since 
E’isjp 8 ~ lOOi^^F 3 , Table 2). This simple argument pre- 
dicts the SGRB flare candidates to be off the LGRB flare 
lag-luminosity relation of a factor > 100 as observed (the > 
inequality accounts for the possible underestimation of the 
real lag due to the fixed and limited energy band- passes bias 
of the previous paragraph). This finding would support the 
presence of non-mesurable (but still non- zero) lags for the 
short burst prompt emission. 

4.7 The flare candidates internal origin 

The above indications point to a link between the properties 
of flare candidates and prompt pulses in SGRBs (for LGRBs 
this was demonstrated by Margutti et al. 2010b). This result 
would naturally arise if both kind of emission share a com- 
mon origin: we refer to this possibility as the internal origin 
interpretation. Theoretical models consistent with this pic- 
ture can be divided into two categories: models requiring a 
late-time GRB central engine activity and models where the 
central engine is not required to be active at the time of the 
flare detection. 

The second class of models includes the magnetic 
re-connection interpretation (Lyutikov & Blandford 2003; 
Giannios 2006): flares would originate from residual, late- 
time magnetic re-connection events within the original out- 
flow (the same ejecta powered the prompt phase) triggered 
by its deceleration due to the sweeping-up of the external 
medium. The deceleration of the original ejecta during the 
afterglow phase causes an increase in the size of causally 
connected regions, thus enabling re-connection of increas- 
ingly larger magnetic structures. The smooth continuum 
would be instead afterglow emission from the shocked ex- 
ternal medium. 

Alternatively, flares and prompt pulses would automat- 
ically share a set of observational properties if they both 
directly arise from the GRB central engine activity (first 
class of models above). If this is the case, the central engine 
would be active on much longer time-scales than previously 
thought (see e.g. Perna et al. 2006; Rosswog 2007; Lee et al. 
2009); at the same time, the similarity of LGRB flares and 
SGRB flare candidates in the L p k/ (L) vs. fpk/Tbo plane as 
well as in the A F/F vs. fpk/Tho space would point to a 
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similar late time evolution of long and short GRB central 
engines. 

It is not possible to observationally discriminate be- 
tween the two scenarios using the present set of data. Careful 
modelling is required (Margutti et al. in preparation). 


5 CONCLUSIONS 

This work presents the first comprehensive analysis of flare 
candidates in SGRBs and provides a comparison to the ob- 
servational properties of X-ray flares in LGRBs with the aim 
to draw an observational picture of SGRB X-ray variability 
any theoretical model is required to explain. 

Our analysis shows that the origin of the SGRB X-ray 
light-curve variability is independent from the large-scale 
host galaxy environment and is not correlated to the de- 
tected afterglow life-time. We find that flare candidates ap- 
pear in different types of SGRB host galaxy environments 
and show no clear correlation with the X-ray afterglow life- 
time; flare candidates are detected both in SGRBs with a 
bright extended emission (EE) in the soft 7-rays and in 
SGRBs which do not show such component (Sec. 4.2). We 
cannot exclude that flare candidates originate on top of faint 
(and undetected in the 7-rays) EEs. In particular, SGRB 
flare candidates are consistent with being drawn from the 
LGRB flare population when considering: 

1. The flare to prompt energy ration £iao,F/Eiao,p (Fig. 
5, Sec. 4.5); 

2. The relative variability time scale A t/t < 1 (Fig.. 7, 
Sec. 4.3): 

3. The rest-frame flare width evolution with time w(t p k) 
(Fig. 2 , Sec. 4.3). 

4. The hard-to-soft trend of the emitted radiation (see 
e.g. Fig. 1). 

The main parameter driving the duration of the episodes 
of variability is the elapsed time horn the explosion, with 
very limited dependence on the different progenitors, envi- 
ronments, life-times, prompt variability time scales and en- 
ergy budgets. The origin of the flare w(f P k) relation must 
arise from what is in common for the long and short burst 
models. From another perspective this result implies that for 
tpk > 100 s the flare duration is likely to retain no memory 
of the variability time-scales of the original prompt mecha- 
nism at work. This would explain why the flare to prompt 
pulse width ratio is different for long and short GRBs. 

On the contrary, SGRB flare candidates significantly 
differ from the standard X-ray flare emission observed in 
LGRBs at the same f P k/(l 4- z) in terms of: 

5. Peak luminosity L^ B ~ O.OlLjgfP (Fig. 3, Fig. 4, 
Sec. 4.5); 

6. Isotropic energy Ef s ® f B ~ O.OlEigo p B (Fig. 5, Sec. 
4.5); 

7. Flare to prompt luminosity ratio Lpk^/Lpk.p- Flare 
candidates in SGRBs are ~ 100 times dimmer than in 
LGRB; 

8. Relative variability flux A F/F (Fig: 7): we find 
A F/F rsj 1 for all SGRB flare candidates (Sec. 4.4); 

9. Lag-luminosity relation: like SGRB prompt pulses, 
flare candidates show shorter lags than expected from the 
lag-luminosity relation of LGRB flares (Fig. 6, Sec. 4.6). 


However and more importantly, the differences listed at 
points 5., 6., 7. and 8. above disappear once the different 
time scale of evolution of the long and short GRB central 
engine as well as the different energy scaling of the two sys- 
tems are properly accounted for (Fig. 8, Sec. 4.5.1 and 4.5.2). 
This finding provides a connection between the properties of 
the detected SGRB X-ray light-curve variability and LGRB 
flares, suggesting a common, internal origin. As a result, 
we conclude that similar dissipation and/or emission mech- 
anisms are likely to be responsible for the prompt and flare 
emission in long and short GRBs, with SGRBs being less en- 
ergetic albeit faster evolving versions of the long category. 
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Figure B4. Same as Fig. B1 for GRB 090607. 
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Table Al- Best fitting parameters of SGRB flare candidates. From left to right: GRB, redshiffc; start time ( t 3 ) y amplitude (A) and shape 
parameters (ri,^) of the best fitting Norris et al. (2005) profile; width evaluated between 1/e intensity points; peak time (£ pk ); relative 
variability time-scale (A t/t = u>/£ p k); relative variability flux A F/F: the value F is computed from the best fit of the continuous emission 
underlying the flare candidates (black solid line of Fig. 1); isotropic equivalent peak luminosity (L pea k) and energy (J£jso) computed in 
the 0.3 — lOkeV observer frame energy band. 


GRB 

z 

ta 

(0 

A 

(c/s) 

ri 

M 

V2 

(») 

w 

(») 

£pk 

(s) 

At/t 

AF/F 

L|P B 

(10 47 ergs -1 ) 

bL°1 e 

(10 48 erg) 

050724 

0.258 

230.0 

12.29 ± 0.88 

12.40 ±3.00 

35.50 ±4.20 

65.11 ±4.64 

250.98 ±1.80 

0.26 ± 0.03 

1.72 ±0.27 

1.49 ±0.34 

7.16 ± 1.98 

0512X0 

- 

107.0 

5.70 ±1.60 

430.00 ± 220.00 

1.77 ±0.87 

14.09 ±3.52 

134.59 ±1.43 

0.10 ±0.03 

1.10 ±0.37 


- 

051227 

- 

105.0 

12.20 ±2.10 

8.90 ± 3.70 

5.90 ± 1.70 

14.35 ± 2.36 

112.25 ±0.88 

0.13 ± 0.02 

1.62 ±0.52 



051227 

- 

122.0 

6.00 ± 0.00 

10.00 ±0.00 

17.40 ± 2.50 

34.94 ±4.13 

135.19 ±0.96 

0.26 ±0.03 

0.80 ± 0.21 



060313 

- 

105.0 

2.80 ±1.20 

19.00 ± 23.00 

4.70 ±5.30 

14.13 ± 8.87 

94.45 ± 2.21 

0.15 ± 0.09 

0.69 ± 0.32 


- 

060313 

- 

150.0 

4.20 ±1.70 

111.00 ±80.00 

10.30 ± 7.10 

38.72 ± 14;65 

183.81 ±4.50 

0.21 ± 0.08 

1.03 ±0.46 



070724A 

0.457 

75.0 

7.80 ±1.70 

25.00 ± 23.00 

7.80 ± 9.00 

22.28 ± 15.76 

88.96 ±2.68 

0.25 ± 0.18 

0.93 ± 0.35 

3.64 ±1.18 

5.07 ±1.64 

070724A 

0.457 

90.0 

12.20 ± 3.90 

42.00 ± 21.00 

5.50 ± 2.10 

19.10 ± 3.60 

105.20 ±1.26 

0.18 ±0.03 

1.46 ±0.65 

5.68 ±2.32 

6.72 ± 2.75 

070724A 

0.457 

150.0 

3.79 ± 0.77 

27.00 ± 17.00 

43.00 ± 18.00 

87.80 ±21.55 

184.07 ±6.34 

0.48 ±0.12 

0.45 ± 0.17 

1.73 ±0.74 

9.73 ±4.16 

071227 

0.383 

150.0 

4.36 ± 0.94 

80.00 ±41.00 

7.30 ± 3. DO 

27.55 ± 5.74 

174.17 ±2.36 

0.16 ±0.03 

0.80 ± 0.29 

1.13 ±0.44 

2.03 ±0.79 

090607 

- 

89.0 

7.60 ±1.10 

53.00 ±31.00 

14.00 ± 12.00 

41 .49 ±22.22 

116.24 ±4.°4 

0.36 ±0.19 

1.54 ±0.45 



090607 

- 

122.0 

3.30 ± 2.90 

18.00 ± 17.00 

39.00 ± 16.00 

75.19 ±25.52 

148.49 ± 12.73 

0.51 ± 0.18 

0.67 ±0.61 



1001m 

0.920 

130.0 

4.90 ±0.73 

2.00 ± 1.80 

30.00 ±0.00 

42.77 ±4.98 

137.75 ±3.55 

0.31 ± 0.04 

0.80 ± 0.56 

9.97 ± 3.48 

21.43 ± 7.48 

100117A 

0.920 

164.0 

7.03 ± 0.94 

6.00 ±0.00 

42.90 ±6.80 

67.78 ±9.13 

180.04 ±1.27 

0.38 ±0.05 

1.14 ±0.79 

14.09 ± 5.60 

47.45 ± 18.86 

100117A 

0.920 

200.0 

4.30 ±1.30 

101.00 ± 77.00 

5.10 ±3.50 

22.11 ±7.86 

222.70 ±2.57 

0.10 ±0.04 

0.70 ± 0.52 

8.49 ±4.57 

8.78 ± 4.73 



Figure El. 0,3-10 keV count-rate light-curve of GRB 050724 and GRB 051210. Black solid line: continuous X-ray emission underlying 
the flare candidates computed as described in Section 2.2; dashed lines: best- fitting flare candidate emission; red solid line: best estimate 
of the total emission. The vertical dot-dashed lines mark the flare candidate onset times. Inset: Complete Swift-XKT light-curve. The 
yellow filled area marks the time window for the computation of the CCF lag. Middle panels: hardness ratio (HR) evolution with time; 
the HR is computed between 1.5-10 keV (hard band) and 0.3-1. 5 keV (soft band). Lower panels: Spectral photon index evolution with 
time as computed by Evans et al. (2010). 
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Figure B2. Same as Fig. B1 for GRB 051227 and GRB 060313. 
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Figure B3. Same as Fig. B1 for GRB070724A and GRB 071227. 


